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Abstract 

We introduce stochastic models of chemotaxis generahzing the deterministic Keher- 
Segel model. These models include fluctuations which are important in systems with 
small particle numbers or close to a critical point. Following Dean's approach, we derive 
the exact kinetic equation satisfied by the density distribution of cells. In the mean field 
limit where statistical correlations between cells are neglected, we recover the Keller-Segel 
model governing the smooth density field. We also consider hydrodynamic and kinetic 
models of chemotaxis that take into account the inertia of the particles and lead to a 
delay in the adjustment of the velocity of cells with the chemotactic gradient. We make 
the connection with the Cattaneo model of chemotaxis and the telegraph equation. 

1 Introduction 

In biology, many organisms (bacteria, amoebae, cells,...) or social insects (like ants, swarms,...) 
interact through the process of chemotaxis [H [21 [3]. Chemotaxis is a long-range interaction 
that accounts for the orientation of individuals along chemical signals that they produce them- 
selves. Famous examples of biological species experiencing chemotaxis are the slime mold 
amoebae Dictyostelium discoideum, the flagellated bacteria Salmonella typhimurium and Es- 
cherichia coli, the human endothelial cells etc. When the interaction is attractive, chemo- 
taxis is responsible for the self-organization of the system into coherent structures such as 
peaks, clusters, aggregates, fruiting bodies, periodic patterns, spirals, rings, spots, honey- 
comb patterns, stripes or even filaments. This spontaneous organization has been observed 
in several experiments [H El El [71 [HI El [IDl [HI [121 [131 [HI [13 HB] and numerical simulations 
[ni[l8l[T9l[2Ql[2Il[22l[23l[2a[2a[26l[2Il[2Hll29l[3a Chemotactic attraction 

is therefore a leading mechanism to account for the morphogenesis and self-organization of 
biological systems. For example, it has been advocated to explain aggregation patterns in 
bacteria, tissue organization during embryonic growth, cell guidance, fish skin pigmentation 
patterning, angiogenesis in tumour progression and wound healing, formation of plaques in 
Alzheimer's disease, dynamics of blood vessel formation etc [211 [M]- It is fascinating to realize 
that the self-organization of chemotactic species in biology shares some analogies with the self- 
organization of galaxies in astrophysics and large-scale vortices (like Jupiter's great red spot) in 



1 



two-dimensional turbulence A first successful model of chemotactic aggregation is provided 
by the Keller-Segel (KS) model [H] introduced in 1970. The standard KS model can be written 
as 

= V-p.Vp-xpVc), (1) 



— = D^Ac- kc + hp. (2) 

It consists in two coupled differential equations that govern the evolution of the density of cells 
(or other biological entities) p(r, t) and the evolution of the secreted chemical c{r,t). The first 
equation ([Tj) is a drift-diffusion equation. The cells diffuse with a diffusion coefficient and 
they also move in a direction of a gradient of the chemical (chemotactic drift). The chemotactic 
sensitivity x is a measure of the strength of the influence of the chemical gradient on the flow 
of cells. The coefficient x can be positive or negative. In the first case (chemoattraction), the 
particles climb the chemical gradient and form clusters. In the second case (chemorepulsion), 
they descend the chemical gradient and repell each other. In that case, the chemical acts like a 
poison. The second equation ([2]) in the KS model is a reaction-diffusion equation. The chemical 
is produced by the bacteria with a rate h and is degraded with a rate k. It also diffuses with 
a diffusion coefficient Dc- When chemotactic attraction prevails over diffusion, the KS model 
describes a chemotactic collapse leading to aggregates or Dirac peaks. There is a vast literature 
on this subject. We refer to Perthame [12] for numerous references in applied mathematics and 
to Chavanis [13] for additional references in physics. 

The first equation of the KS model can be interpeted as a mean-field Smoluchowski equation 
describing a system of Brownian particles in interaction. On the other hand, in the limit of 
large diffusivity of the chemical, we can make a quasi-stationary approximation dc/dt ~ in 
the second equation and obtain the screened Poisson equation. We are led therefore to the 
simplified Keller-Segel model 

^ = V-iD,Vp-xpVc), (3) 



Ac — /cqC = — Ap, 



(4) 



where we have set 



k/Dc and A = h/Dc- In the absence of degradation of the chemical 
{ko = 0), the field equation (j4]) reduces to the Poisson equation Ac = — Ap (see [44] and 
Appendix C of [32] for a precise justification of these approximations). In that case, the Keller- 
Segel (KS) model becomes isomorphic to the Smoluchowski-Poisson (SP) system 



dt 



(5) 



A* = SjGp, (6) 

^Tliese analogies arc intrinsically due to the long-range attractive nature of the interaction. In particular, self- 
gravitating systems, 2D vortices and chemotactic species interact through a field produced by the distribution of 
particles via a Poisson equation (or its generalizations). Furthermore, the process of self-organization is described 
by relatively similar relaxation equations corresponding to nonlinear mean field Fokkcr-Planck equations |35j . 
Therefore, self-gravitating systems, 2D vortices and chemotactic species share many analogies despite their 
very different physical nature. These striking analogies have been emphasized by the author in several papers 

[aelEilEHllMlliniES]. 
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describing a system of overdamped self-gravitating Brownian particles in the mean field approx- 
imation [I5l[25l[28l[30lll6l[nillllll3lll8lll9]. We have the correspondances: = kBT/^m, 
X = 1/^, c = — $, A = SdG. In particular, the concentration of the secreted chemical 
c(r, t) = — $(r,t) in biology plays the role of the gravitational potential (with the opposite sign) 
in astrophysics^. More generally, when we consider a system of Brownian particles interacting 
via an arbitrary binary potential u{r — r') and make a mean- field approximation 
we obtain the mean-field Smoluchowski equation 



dp 
di 



V 



m 



Vp + pV$ j 



(7) 



^r,t) = J p{r',t)u{r-r')dr'. (8) 

The main difference between models ([I])-© and ©-([H]) comes from the equation for the field 
c{r,t) or $(r, t). Equation ([2]) is non-markovian since the concentration of the chemical c{r,t) 
at time t depends on the concentration of the bacteria and of the chemical at earlier times. 
By contrast, Eq. ([8]) is markovian since the potential $(r,t) is assumed to be instantaneously 
produced by the distribution of particles. 

It is important to note that the Keller-Segel model is a mean field model which ignores 
fluctuations. This implicitly assumes that the number of cells — > -|-oo and that we are far 
from a critical point [56]. Now, in biology, the number of particles in the system can be relatively 
small. Furthermore, from the statistical physics viewpoint, it is natural to investigate the role 
of fluctuations during chemotaxis. In order to go beyond the mean fleld approximation, some 
authors [ITJ [STJ [58], [32] have proposed to return to a corpuscular description of the dynamics 
and to describe the motion of the particles (chemotactic species or "active" walkers) by N 
coupled stochastic Langevin equations of the form 

^ = xVQ(r,(t),t) + v^Ri(t), (9) 



= D.Acd -kc, + hJ2 <^(r - ri{t)), (10) 

i=l 

where rj(t) denote the positions of the particles, Cd{r,t) is the exact fleld of secreted chemical 
and Hiit) is a white noise satisfying (Ri(t)) = and {Ri^ait)Rj^i3{t')) = 5ij5ap5{t — t') where 
i = 1,...,N refer the the particles and a = l,...,d to the dimensions of space. Note that the 
motion of cells is treated on an individual basis but the chemical signals are treated in the 
continuum limit. This separation of scales appears to be reasonable in most applications. In 
the mean fleld approximation, these stochastic equations lead to the KS model ([I])-([2])|fl. When 
the reaction-diffusion equation ( |T0|) is replaced by a Markovian equation of the form 

N 

AQ-A:2crf = -A^5(r-r,(t)), (11) 

1=1 

^One great achievement of Keller & Segel [H] was to interpret slime mold aggregation as a manifestation of 
a fundamental instability in a uniform distribution of amoebae and acrasin (chemoattractant). As noticed in 
[501 [5T] , this instability is closely related to the Jeans gravitational instability in astrophysics [Slj . 

■^Stevens |57| gives the first rigorous derivation (in the mathematical sense) of the KS model from an inter- 
acting stochastic many-particle system where the interaction between the particles is rescaled in a moderate 
way as the population size N tends to infinity. 
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we obtain a simplified model of chemotaxis that leads to the simplified KS model ©-dlD in 
the mean field approximation. More generally, for Brownian particles interacting via a binary 
potential of interaction u{r — r'), one obtains the stochastic model 



considered in [SHI EDI EB ESI ESI [Ml SZl ES] • In the mean field approximation El ES], 
these equations yield the mean- field Smoluchowski equation (I7])-([HD. 

In systems with weak long-range interactions, the mean field approximation is expected 
to become exact in a proper thermodynamic limit — >■ -|-oo such that the strength of the 
potential scales like while the volume V remains of order unity [51] • In the context of 
chemotaxis, the differences between mean field and non mean field models have been discussed 
by Grima [3lj who showed situations where the mean field approximation fails to predict the 
width of the aggregate sizes. In particular, the disagreement is very severe close to the critical 
point where we know that mean field approximation breaks down in general [63]. This is 
because the fluctuations become very important so that it is not possible to neglect the two- 
body correlation function anymore [56]. On the other hand, the mean field approximation 
assumes that the number of particles ^ 1. In stellar systems and plasmas, this is always the 
case. However, for biological systems, the number of interacting bacteria or cells is frequently 
less than a few thousands so that finite N effects and statistical fluctuations are important. 
In view of these remarks, it is highly desirable to obtain stochastic kinetic equations that take 
into account fluctuations and that go beyond the deterministic mean field Keller-Segel model. 
Such equations are discussed in the present paper. In the first part of the paper (Sec. [2]), 
following Dean's approach [HD], we derive the exact kinetic equation satisfied by the density 
distribution of chemotactic species. This equation takes into account stochastic fluctuations 
and memory effects present in the field equation for the secreted chemical. If we average over 
the noise, we recover the hierarchy of kinetic equations discussed by Newman & Grima [58] . 
If we make a mean- field approximation, we recover the Keller-Segel model [H]. Therefore, 
this exact stochastic kinetic equation generalizes several models introduced in the chemotactic 
literature. We also propose a simplified kinetic equation for a coarse-grained density field 
(instead of a sum of (5-functions) keeping track of fluctuations. This equation (l3Tl) - (!32l) could 
be of practical interest in chemotaxis. In the second part of the paper (Sees. [3]toE]), we note 
that the Keller-Segel model is a parabolic model which neglects the inertia of the particles 
and which assumes an instantaneous adjustment of the velocity with the chemotactic gradient. 
We consider hyperbolic models that generalize this parabolic model. We first consider the 
Cattaneo model of chemotaxis [26] which consists in introducing a delay in the establishment 
of the current (Sec. [3]). Then, we consider hydrodynamic models including a friction force (Sec. 
H]). Using a semi-linear approximation, we show that the Cattaneo model can be recovered from 
these hydrodynamic equations [SS]- In Sec. El we generalize these models so as to take into 
account fluctuations. This leads to stochastic hyperbolic models of chemotaxis which generalize 
the ordinary deterministic parabolic Keller-Segel model. Finally, in Sec. El we develop a kinetic 
theory of chemotactic species in phase space taking into account the inertia of the particles and 
the discrete nature of the system. We derive stochastic kinetic equations that should improve 
the description of the cells' motion. The link with the parabolic and hyperbolic models is also 
discussed. 




(12) 



N 




(13) 



i=l 
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This paper adapts the results of [SB] to the context of chemotaxis with complements and 
amplification. Although these different stochastic equations (in particular the parabolic ones) 
are well-known in statistical physics [591 EOl ED [62], [56] their application to the context of 
chemotaxis, proposed in [56], is new and is an important contribution of the present paper. 



2 The stochastic Keller-Segel model 



In this section, we introduce a stochastic model of chemotaxis, generalizing the Keller-Segel 
model, by taking into account fluctuations. Let us first derive the exact kinetic equation satisfied 
by the density distribution of cells whose dynamics is described by the coupled stochastic 
Langevin equations (p|)- (fTO|l . We follow Dean's approach [60]. The exact density field, expressed 
in terms of 5-functions, can be written 



TV 



TV 



(14) 



1=1 



i=l 



For any function -F(r), we have F(rj(t)) = J pj(r, t)F(r)(ir. Now, using Ito's calculus [51], one 
has 



dt 



p,(r, t) xVF(r) ■ VQ(r, t) + y/2D,VF{r) ■ R,{t) + D,AF{r) 



dr. 



(15) 



Integrating by parts, we obtain 
dF{r,) 



dt 



F{r) -xV-(p.(r,t)VQ(r,t))- v^V-(p,(r,t)Ri(t)) + D,Ap,(r,t) rfr. (16) 



Then, using dF{ri)/dt = f dtPi{r,t)F{r)dr and comparing with Eq. (fT6l) . we get (using the 
fact that F is an arbitrary function) 



dpi 
dt 



-XV ■ (pi(r, t)VQ(r, t)) - ^2D,W ■ {pi{r, t)Ri{t)) + D,Api{r, t). 
Summing this relation over the i, we finally obtain 

dpd, 



(17) 



TV 



dt 



[r, t) = D,Apd{r, t) - xV ■ (pd(r, t) VQ(r, t)) - v^$^ V ■ (pi(r, t)K,{t)). (18) 



i=l 



Now, the last term can be rewritten 1601: 



N 



5^V- (p.(r,t)R,(t)) = V - (py'(r,t)R(r,t)), 



(19) 



i=l 



where R(r, t) is a Gaussian random field such that (R(r, t)) = and {Ra{jc,t)Rp{Y' ,t')) = 
5a/3'^(r — r')(5(t — t'). Therefore, the system of equations satisfied by the exact density field 
expressed in terms of (5-functions is 



dpd 
dt 



r, t) = D.Apdir, t) - xV ■ (pd(r, t) VQ(r, t)) + V ■ ( V2Z},pd(r, t)R(r, t) 



(20) 



dcd 
dt 



= D,Acd{r, t) - kcd{r, t) + hp^iv, t) . (21) 
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The first and third terms in the r.h.s. of Eq. fl20p correspond to a pure Brownian motion 
and the second term takes into account chemotaxis, i.e the attraction or repulsion of the cells 
by the chemical. As noted by Dean [60], the noise in Eq. (120|) appears not additively but 
multiplicatively. 

Integrating Eq. fllOp . the concentration of the chemical can be expressed in terms of the 
cell paths as 



Cd 



p rt N 

'r, t) = h dr' dt'Giv -r',t- t') V 5{r' - v,{t')) 
J i=i 



(22) 



where the Green function for the chemical diffusion equation is given by 



^(r,^) = {ATiD.ty/^^v -TTTT - kt 

The gradient of the concentration field is 

Vcrf(r,t) = /i j dr' dt'VG{r -r',t-t')pd{r',t' 

Substituting Eq. 1^ in Eq. ([20D, we obtain 
dpd 



(23) 



(24) 



dt 



{r,t)=D,Apdir,t)-xhV 



Pd{r,t)V j dr' dt'Gir -r',t-t')pa{r',t' 
+V- (V2Z},pd(r,t)R(r,t) 



(25) 



If we average over the noise and introduce the smooth density p{r,t) = {pd{r,t)), we recover 
Eq. (9) of Newman & Grima [58] : 

^{r,t)=D,Ap{r,t)-xhV- j dr' j\t'[VG{r - r' ,t - f)]{pa{r,t)p,(r' ,t')) . (26) 

If we make a mean field approximation {pdi^, t)p^(r', t')) ~ p(r, t)p(r', t') in Eq. (126|) . we recover 
the Keller-Segel model 



with 



|^(r, t) = D*Ap(r, t) - xV ■ (p(r, t) Vc(r, t)), 



c{r,t) = h J dr' dt'G{r-r',t-t')p{r',t'). 



(27) 



(28) 



Given the definition of the Green function G, the smooth concentration c(r, t) is solution of the 
reaction-diffusion equation 



Oc 

— = DcAc — kc + hp. 
dt ^ 



(29) 



We also note, for future reference, that the steady solutions of the KS model fl27j) correspond 
to a mean field Boltzmann-like distribution 



P 



(30) 
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where T^ff = D^/x is an effective temperature given by an Einstein relation. 

Grima [M] has shown that the mean field approximation may lead to wrong results if we 
are close to a critical point or if the number of particles is not large enough. Therefore, it 
may be useful to have a more general model than the Keller-Segel model ([I])-© which keeps 
track of fluctuations. Equations fl20p - (pT]) are exact and contain the same information as the 
A^-body stochastic Langevin equations (P)- (niII) . They are not very useful for practical purposes 
since they govern the evolution of a density field which is expressed as a sum of 5-functions. 
It is easier to directly solve the equivalent A^-body stochastic Langevin equations fl9|l- f|T0|) . 
However, using phenomenological arguments like those described in [62l [56] , we can consider a 
spatio-temporal coarse-grained distribution p(r, t) which smoothes out the exact density field 
Pd{T',t) while keeping track of fluctuations. We also assume that the spatio-temporal window 
is sufficiently small so that we can make the approximation p^^''(r, r', t) ~ p(r, t)p(r', t). In that 
case, we obtain the stochastic Keller-Segel model for the coarse-grained distribution 

|^(r, t) = D,Ap(r, t)-xV- (p(r, t) Vc(r, t)) + V- (v/2Z},p(r, t)R(r, t)) , (31) 



dc 

— = DcAc — kc + hp, 
dt ^' 



(32) 



generalizing the deterministic Keller-Segel model ([I])-©. This equation is one of the most 
important result of this paper. As shown in Appendix B of [56], the form of the noise term 
in Eq. (l3T|l can be obtained from the general theory of fluctuations developed in Landau & 
Lifshitz [65]. This provides another, direct, justification of the stochastic Eq. (!3T]) . As shown 
in [56], the mean field approximation breaks down close to a critical point because the two- 
body correlation function diverges. In that case, it may be more relevant to use the stochastic 
Keller-Segel model f l31l) - fl32l) including fluctuations instead of the deterministic Keller-Segel 
model (HI)-®. 

It is also very important to take into account fluctuations when the system can be found in 
several metastable states. If we introduce the coarse-grained free energy functional 



X 



J plnpdr + ^ j [DcC^cY + kc^] dr - J p cdr, 



(33) 



we can write the stochastic equation (13T!) in the form 



dt 



Xp(r,t)V-^^^-^- 



6p 



V- 



2D,p(r,t)R(r,t) 



(34) 



This equation can be viewed as a Langevin equation for the field p(r, t). The evolution of the 
probability of the density distribution iy[p, t] is governed by a Fokker-Planck equation of the 
form 



dW[p,t] 
dt 



6p{r, t) 



V-p(r,t)V 



5 
-Tp 



6p 



W[p,t] \dr. 



(35) 



At equilibrium, we have W[p] oc e 



'Z^'^'" with Fc.g.[p] = ^ Jplnpdr 



Up cc^r 



(we have substituted Eq. fl32|) with dc/dt = in Eq. fl33|) ). For A^ +oo, the equilibrium 
distribution W[p] is strongly peaked around the global minimum of Fc,g,[p] at fixed mass M = 
J J) dr. However, the system can remain trapped in a metastable state (local minimum of 
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F(.,g\p\) for a very long time which becomes infinite at the thermodynamic hmit +00. Let 
us be more precise. If we ignore the noise term, Eq. (!34l) reduces to 



dt 



(36) 



5-p 

which is the deterministic Keller-Segel model ( l27l) . This equation satisfies an H-theorem 

F = - — (D^Vp - xP^cf dr-^ (DcAc -kc + hp? dr < 0, (37) 
J XP hj 

with F = iff the distribution is given by Eq. ( l30l) . Therefore, a steady state is stable iff it is a 
(local) minimum of free energy at fixed mass. Assuming that the free energy is bounded from 
below, we know from Lyapunov's direct method that the system will relax towards a steady 
state that is a minimum (global or local) of the free energy functional Fc_g,[p] at fixed mass 
(maxima or saddle points of free energy are linearly dynamically unstable with respect to mean 
field Fokker-Planck equations [35]). If the free energy admits several local minima, the selection 
of the steady state will depend on a notion of basin of attraction. Without noise, the system 
remains on a minimum of free energy forever. Now, in the presence of noise, the fiuctuations 
can induce dynamical phase transitions from one minimum to the other. We should therefore 
see the system "jump" between different states. Thus, accounting correctly for fluctuations 
is very important when there exists metastable states. The probability of transition scales 
as e~^^^'^^fi' where AF is the barrier of free energy between two minima. Therefore, on an 
inflnite time, the system will explore all the minima and will spend most time in the global 
minimum. This will be the case only if is not too large. Indeed, for systems with long-range 
interactions, the barrier of free energy AF scales like A^ so that the probability of escape from 
a local minimum is very small and behaves like . Therefore, even if the global minimum is 
in principle the most probable state, metastable states are highly robust in practice since their 
lifetime scales like . They are thus fully relevant for A^ ^ 1: metastable states are in practice 
"stable states". These interesting features (basin of attraction, dynamical phase transitions, 
metastability,...) would be interesting to study in more detail in the case of chemotaxis. The 
study of the stochastic Keller-Segel model will be considered in future publications. 

3 The Cattaneo model of chemotaxis 

The general Keller-Segel (GKS) model [H] can be written as 

= V ■ (D2(P, c) Vp - Di(p, c) Vc), (38) 



— = D,Ac-k{c)c + h{c)p, (39) 

where Di = Di{p,c) and D2 = D2{p,c) can both depend on the concentration of the cells 
and of the chemical. This takes into account microscopic constraints, like close-packing effects, 
that can hinder the movement of cells and lead to nonlinear diffusion and nonlinear mobility. 
The GKS model 0381) - fl5^ can be viewed as a nonlinear mean field Fokker-Planck equation 
associated with a notion of effective generalized thermodynamics [35]. The flrst equation can 
be written in the form of a continuity equation dtp = —V ■ J with a current 

J = {D^ip, c)Vp - D,{p, c)Vc) . (40) 
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It is important to note that the GKS model is a parabolic model like the usual heat diffusion 
equation. Like for the Fourier law of heat conduction, it is assumed that the current J is 
instantaneously equal to the right hand side of Eq. (l40l) . that we shall call the "chemotactic 
gradient" for future reference. In the context of heat conduction, Cattaneo [66] has proposed 
a modification of Fourier's law in order to describe heat propagation with finite speed. In the 
context of chemotaxis, Dolak & Hillen [25] have introduced a Cattaneo model for chemosensitive 
movement. They assume that the current is not instantaneously equal to the chemotactic 
gradient but relaxes to it with a time constant 1/r. Then, the corresponding Cattaneo model 
for chemosensitive movement reads 



Taking the time derivative of Eq. (HTj) and using Eq. (l42l) . we obtain the hyperbolic model 



This equation, which is second order in time, is analogous to the telegraph equation which 
generalizes the diffusion equation by introducing memory effects. For r = 0, we recover the 
GKS model (138D-([39D particular case. 

4 Hydrodynamic models of chemotaxis 

The parabolic Keller-Segel model [H] is able to reproduce the formation of clusters (clumps) 
resulting from chemotactic collapse. This can explain experiments on bacteria like Escherichia 
coli or amoebae like Dictyostelium discoideum exhibiting pointwise concentrations [H [HI [9l [HI 
[SI [6l [71 [8] • Recently, several experiments with human endothelial cells have shown the formation 
of networks that can be interpreted as the initiation of a vasculature [TOl [131 [23 [ISl [15]. Cells 
randomly spread on a gel matrix autonomously organize to form a continuous multicellular 
network which can be described as a collection of nodes connected by chords [27]. This process 
takes place during the early stages of vasculogenesis in embryo development. These filaments 
are observed in the experiments of capillary blood vessel formation. These structures cannot be 
explained by the Keller-Segel parabolic model which generically leads to pointwise blow-up 0. 
In order to account for these filaments, hyperbolic models of chemotaxis have been introduced 
[271 [SZl [Ml [23 [SDl [221 [SB] • They have the form of damped hydrodynamic equations taking 

^In fact, some Keller-Segel models including cell kinetics can, under certain conditions, give rise to network- 
like patterns (see Fig. 12c of [24]). 

^The type of hydrodynamic equations (I44ll - (l45p including a long-range mean field interaction, a density 
dependent pressure and a friction force were introduced in Chavanis [69j (see also [36]) for Langevin particles in 
interaction and called the damped Euler equations. Their application to chemotaxis and gravity was mentioned. 
These equations can be derived from kinetic equations (nonlinear mean-field Fokker-Planck equations) by using a 
local thermodynamic equilibrium condition (L.T.E.) to close the hierarchy of hydrodynamic moments [691147115^ . 
However, they remain heuristic because the L.T.E. approximation is not rigorously justified. By contrast, in 
the strong friction limit ^ -\-oo, we can rigorously derive the GKS model (|55)) - ([55)) . also called the generalized 
Smoluchowski equation, by using a Chapman-Enskog expansion |70| or a method of moments [711 1321 135] . The 
model considered by Gamba et al. [27] (see also 67j2^2M) corresponds to ^ = in Eq. (|45)) . It can be derived 
in an asymptotic limit of kinetic equations of a different type (see [29] and Appendix D of [32]). In more recent 
papers [68], the aforementioned authors have also included a friction force in their model. 




(41) 




(42) 




(43) 



9 



into account inertial effects : 

1^ + V ■ (pu) = 0, (44) 

d 

— (pu) + V (pu ® u) = -D2 (p, c) Vp + Di (p, c) Vc - epu. (45) 

Considering the momentum equation P5|) . the inertial term (l.h.s.) models cells directional 
persistence, i.e. the natural tendency of a particle to continue in a given direction in the absence 
of any interaction. When D2{p, c) depends only on the density, the first term on the r.h.s. can 
be interpreted as a barotropic pressure force — Vj9(p) (see [35] for different examples of equations 
of state). The pressure law is expected to be linear for low densities and to increase rapidly 
above a certain threshold ~ ctq in order to describe the fact that the cells do not interpenetrate. 
For example, in [691 HOI [32] we proposed to take p(p) = — (ToTe//ln(l — p/o"o) which returns 
the "isothermal" equation of state p = pTeff for dilute systems p <^ ctq where the motion of 
an individual cell is not impeded by the other cells [25] , and which diverges when the cells are 
compressed towards the maximum density p ^ o"o. Another possible equation of state is the 
polytropic one p{p) = Kp"' [72l [32] taking into account anomalous transport (normal transport 
corresponds to the isothermal case 7 = 1). The chemotactic response -Di(p, c) of the bacterium 
to the chemical gradient (second term in the r.h.s. of Eq. psj) ) can also depend on c and p so as 
to take into account anomalous reactivity (the normal case corresponds to -Di(p, c) = p but the 
form Di{p, c) = p{l — p/ctq) has also been considered to take into account volume filling effects 
[211 [691 [101 [35]). Finally, the last term in the r.h.s. of Eq. (145|) is a friction force that measures 
the importance of inertial effects. It parametrizes the tendency of the organisms to continue 
in a given direction. In this inertial model, the velocity of a particle takes a finite time to 
get aligned with the chemotactic gradient while in the Keller-Segel model, this alignement is 
assumed to be instantaneous (see below). The "delay" in the alignement of the velocity with 
the chemotactic gradient is similar to the idea that is at the heart of the Cattaneo model in 
Sec. [1 

If we neglect the friction force (^ = 0) we recover the model introduced by Gamba et al. [27] . 
Alternatively, if we neglect the inertial term (l.h.s.) in Eq. fH5l) and substitute the resulting 
expression [521 [SDl [22] : 

pu = -i (Z}2(p, c)Vp - Di(p, c)Vc) , (46) 

in Eq. fl44p . we recover the GKS model. This is valid in a strong friction limit S, +00. We 
can also obtain a more general model taking into account some memory effects. If we neglect 
only the nonlinear term V(pu ® u) in Eq. (145|) . we obtain 

^ (pu) = -1^2 (p, c) Vp + (p, c) Vc - epu, (47) 

which is equivalent to the Cattaneo model fl42p with r = 1/^. Taking the time derivative of Eq. 
OH]) and substituting Eq. P71) in the resulting expression, we obtain a simplified hyperbolic 
model keeping track of memory effects 

|f + = V ■ {D,{p, c)Vp - Di(p, c)Vc) . (48) 

This provides a new justification (see also [56]) of the Cattaneo model of chemotaxis from the 
damped hydrodynamics equation (l4^ - (H5l) . This can be viewed as a semi-linear hydro dynamic 
model since its derivation assumes that the nonlinear term V(pu ® u) in Eq. (H5l) can be 
neglected while the full nonlinearities in the r.h.s. of Eq. ( H5l) are taken into account. 
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5 Stochastic hydrodynamic models of chemotaxis 



In this section, we generalize the previous hydrodynamic equations in order to take into account 
fluctuations. We restrict ourselves to the standard situation where D2 = C,D^ and Di = p. The 
stochastic damped Euler equations generalizing Eqs. (I44l)-( l45|) can be written 

^ + V ■ (pu) = 0, (49) 



d 

— (pu) + V(pu ® u) = -^D,Vp + pVc - ^pu - y/2D,ep R(r, t). (50) 

As shown in Appendix B of [56], the form of the noise in these equations can be obtained by 
applying the general theory of fluctuations developed by Landau & Lifshitz [65]. If we neglect 
the inertial term (l.h.s.) in Eq. (!50|) and substitute the resulting expression 



pu = -{D.Vp - xp^c) - ^2D,p R(r, t), (51) 

where x = 1/^ in Eq. (l49l) . we recover the stochastic Keller-Segel equation (I3T!) . This is valid 
in a strong friction limit ^ — +00 with ~ 1. As in Sec. HJ we can obtain a more general 
model taking into account some memory effects. Indeed, if we neglect only the nonlinear term 
V(pu ® u) in Eq. fl50l) . we find 

X^(pu) = -D,Vp + xpVc - pu - ^/2D7p R(r, t). (52) 

Taking the time derivative of Eq. (149!) and substituting Eq. (152|) in the resulting expression, 
we obtain the stochastic Cattaneo model of chemotaxis 

+ 1^ = V ■ {D.Vp - XP^c) + V ■ {,/2D7pK). (53) 

6 Stochastic kinetic models of chemotaxis 

In order to take into account fluctuations in a rigorous way, we must start from a microscopic 
description of the dynamics of the chemotactic species. In Sec. [21 we have considered an 
overdamped dynamics. However, according to recent observations in biology (as discussed in 
Sec. H]), it is important to take into account the inertia of the particles. A kinetic model 
of chemotaxis taking into account finite N effects and inertial effects has been proposed in 
Chavanis & Sire [32] • In the simplest case, the motion of the biological entities is described by 
coupled stochastic Langevin equations of the form 



-jf = -evi + Vcd{ri{t),t) + V2DIU{t), 
at 



(55) 



dcd 

at 



N 



(56) 



i=l 
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where ^ is a friction coefficient and D a diffusion coefficient in velocity space. We can intro- 
duce an effective temperature Tg// through the Einstein relation T^ff = D/^ [32l [35]. The 
overdamped stochastic equations (p|)- (fTO|) can be recovered in a strong friction limit ^ — > +cx), 
neglecting the inertial term in Eq. (!55|) . and writing X = 1/^ ^ind D,, = D/C,^. We now pro- 
ceed in deriving the exact kinetic equation satisfied by the distribution function of cells whose 
dynamics is described by the coupled stochastic Langevin equations fl54l) - fl56l) . The exact dis- 
tribution function, expressed in terms of (5-functions, can be written 



N 



N 



/,(r, V, t) = ^ /,(r, V, t) = 5^ 6{r - r,(t))5(v - v,(t)). 



(57) 



i=l 



i=l 



For any function F(r, v), we have F(rj(t), Vi(t)) = J /j(r, v, t)F(r, v)(ir(iv. Now, using Ito's 
calculus, one has 



dF{ri,Yi) 
dt 



/i(r, v,t) 



VrF(r, v) ■ V - eVvF(r, v) ■ v + VvF(r, v) ■ VQ(r, t) 



+V2D\/^F{t, v) ■ Ri{t) + DAvF(r, v) 



drdv. (58) 



Integrating by parts, we obtain 
dF{ri, Vi) 



dt 



F r,v 



df d df- 

V ■ ^(r, V, t) + ■ (/.(r, V, t)v) - VQ(r, t) ■ ^(r, v, t) 



d 



'2D— (Mr, V, t)R,,(t)) + DA^Mr, v, t) 



c/rrfv. (59) 



Then, using dF{ri,\i)/dt = J dtfi{r,\,t)F{r,Y)drd'v and comparing with Eq. fl59l) . we get 



^+v- — + Vcd- — = — ■ + ^/iV 

ot or ov av \ av 

Summing this relation over the i, we finally obtain 



d 



'2D— . (/,R0 . 



(60) 



(61) 



Now, proceeding like in [50], the last term can be rewritten: 

- ■ (/.(r,v,t)R.(t)) = ^ ■ (/]/^(r, v,t)Q(r, v,t)). 



i=l 



(62) 



where Q(r, v, t) is a Gaussian random field such that (Q(r, v, t)) = and {Qa{r, v, t)Qp{jc' , v', t')) 
= 5a/35(r — r')5(v — V)5{t — t'). Therefore, the system of equations satisfied by the exact dis- 
tribution function expressed in terms of (5-functions is 



(63) 



— = DcAcd -kcd + h I fd{r, v, t)dv. 



(64) 



12 



This will be called the stochastic Kramers equation of chemotaxis for the exact distribution 
function. Using Eq. (12^ . it can be written 

^ + V . ^ + A y dv'dV f dtVG(r - r', i - «') A(r', v', i') ■ ^(r. v, i) 

If we average over the noise and introduce the smooth distribution function /(r, v, t) = (/^(r, v, t)), 
we recover Eq. (60) of Chavanis & Sire |32j : 



1^ + V . 1^ + ft . j dr'dV f dt'VG(v -v'.t- t'){U(T, V, «)/j(r', v', i')) 



|_.(z,|h-«/v|. («) 



If we make a mean field approximation (/^(r, v, t)fd{r', v', t')) ~ /(r, v, t)f{r', v', t'), we recover 
Eqs. (66)- (68) of Chavanis & Sire [32]: 

^ + v.^ + Vc-^ = — ■ f^^+e/vl (67) 
dt dr dv dv \ dv / ' 



(9c 
di 



D^Ac-kc + h J /(r,v,t)dv. (6J 



This can be viewed as a mean field Kramers equation of chemotaxis in the same way that the 
Keller-Segel model can be viewed as a Smoluchowski equation of chemotaxis. In fact, the Keller- 
Segel model ([I])-([2]) can be recovered from Eqs. (l67I) - (!68l) in a strong friction hmit ^ — * +oo 
by using a Chapman- Enskog expansion [70] or a method of moments [32]. Let us note, for 
future reference, that the steady solutions of the mean field Kramers equation of chemotaxis 
correspond to a mean field Maxwell-Boltzmann-like distribution 

/ = A'e-'^^^''/^--), (69) 

where (3 = l/T^ff is the inverse effective temperature. If we integrate this distribution over 
the velocitities we recover the distribution fl5Ul) that is the steady solution of the Keller-Segel 
model dH)-®. 

As discussed in the Introduction, the mean field approximation may not always give a 
good description of the dynamics. On the other hand, Eqs. (!63|) -( l64|) for the distribution 
function expressed in terms of (5-functions are exact but they are too complicated for practical 
purposes because they contain exactly the same information as the A^-body stochastic Langevin 
equations (l5^ -( l56l) . Therefore, as in Sec. El we shall introduce a simplified kinetic equation 
for a coarse-grained distribution function /(r, v,t) which smoothes out the exact distribution 
function /^(r, v, t) while keeping track of fluctuations. We propose the simplified stochastic 
model 

dj dj dj d ( Bj - \ d 



dc 

— = D,Ac -kc + h I /(r, V, t)dw. (71) 
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This model takes into account inertial effects and fluctuations so that it should provide a good 
description of the dynamics of chemotactic species. As shown in Appendix B of [56], the form 
of the noise in these equations can be obtained by applying the general theory of fluctuations 
developed by Landau & Lifshitz [65] . 

Let us try to make a connexion with the hydrodynamic equations introduced phenomeno- 
logically in Sec. O Taking the hydrodynamic moments of the stochastic Kramers equation (1701) 
and proceeding as in [32], we obtain 

1^ + V ■ (pu) = 0, (72) 

di^^'''^ ^ d^j^^'^''''^ ^~~d^^ ^d^i " ^^""^ ~ J v/2a7<3«c?v, (73) 

where p{r,t) = J fdw is the density, u(r,t) = (1/p) / /vrfv is the local velocity, w = v — 
u(r, t) is the relative velocity and Pij = J fwiWjdv is the pressure tensor. Defining g(r, t) = 
J i/ZD/Qrfv, it is clear that g is a Gaussian noise and that its correlation function is 

{g,{T,t)g^{r',t')) = 2D J ^ f{r, v, t)/(r', v', t') {Q,{r, v, t)Q,{r', v', t'))d^rdV 

= 2D5ij5{r - v')5{t - t') j /(r, v, t)dv = 2D5ij5{v - r')5(t - t')p(r, t)- (74) 



Therefore, the equation for the momentum f[73l) can be rewritten 

This equation is not closed since the pressure tensor depends on the next order moment of 
the velocity. If, following [32], we make a local thermodynamic equilibrium (L.T.E.) approx- 
imation V, t) ^ {(3 p{v,t)e~^^ to compute the pressure tensor, we find that 
Pij ~ Teffp6ij. In that case, Eqs. (!72|) and (!75l) return the stochastic damped Euler equa- 
tions (l49l) - (!50l) . We recall, however, that there is no rigorous justification for this local ther- 
modynamic equilibrium approximation. Therefore, it does not appear possible to rigorously 
derive the damped hydrodynamic equations (HU]) - (|5U]) from the Kramers equation (ITOl) - !!?!!) 
by a systematic procedure. Alternatively, if we consider the strong friction limit ^ +oo 
for fixed (3, implying D = ^/ (3 ^ +oo, the first term in the r.h.s. of Eq. flTOj) implies that 
/(r, v,t) ^ {P/2TiYl^p{v,t)e-^^'l^ + 0(1/0, u = 0(1/0 and P,,- = T,ffp5,, + 0(1/0 PI- To 
leading order in 1/^, Eq. fl75]) becomes 

pM^-j{T,ffS/p- pVc+ ^/2D'pYi{r,t)y (76) 

Inserting Eq. fl75]) in the continuity equation fl72p and recalling that Tg/j = D/C, = ■C-D* 
and X = 1/C; W6 recover the stochastic Keller-Segel model (13T|) - (132|) . It is therefore possible 
to rigorously derive the stochastic Keller-Segel model (13T]) -( !32|) from the stochastic Kramers 
equation (170l) - (17T!) in the strong friction limit ^ — ^ +oo. 

7 Conclusion 

In this paper, we have derived generalized Keller-Segel models of chemotaxis taking into ac- 
count fluctuations. This leads to stochastic kinetic equations instead of deterministic equations. 
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Fluctuations become important close to a critical point [031 EH EH], so it is valuable to have 
a model of chemotaxis going beyond the mean field approximation and taking into account 
fluctuations. The divergence of the spatial correlation function close to the critical point has 
been analyzed in detail in [56] for Brownian particles interacting through a binary potential. 
These particles are described by a stochastic Smoluchowski equation coupled to the markovian 
fleld equation (jS]). The general methods developed in [SB] can be extended to the stochastic 
Keller-Segel model f l3T|) coupled to the non-Markovian field equation ([2]). Accounting for fiuctu- 
ations is also important when the number of particles N is small and when there exists several 
metastable states. In that case, fluctuations can trigger dynamical phase transitions from one 
state to the other. 

We have also introduced kinetic models of chemotaxis in phase space taking into account 
inertial effects. In the strong friction limit, we recover the Keller-Segel model describing an over- 
damped dynamics. We have discussed the relation between the kinetic equations in phase space 
and the hydrodynamic equations introduced phenomenologically. Finally, we have shown how 
the Cattaneo model of chemotaxis [SB] could be obtained from these hydrodynamic equations. 

This paper and [SH] are the flrst attempts to include fluctuations in the kinetic equations of 
chemotaxis (the main results were given in [5B] and they have been discussed here speciflcally 
with more details and amplification). In view of the importance of the Keller-Segel model in 
biology, the stochastic equations that we propose can have a lot of applications and can open 
the way to many new investigations. Their detailed numerical and analytical study is therefore 
of considerable interest. We hope to come to these problems in future works. 

Note added: Until now, fiuctuations have been ignored by people working on chemotaxis. 
Therefore, Ref. [56] and the present paper are the first attempts to include fiuctuations in 
the Keller-Segel model. However, after submission of these papers, a paper by Tailleur & 
Gates [arXiv: 0803711369] (now pubhshed as Phys. Rev. Lett. 100, 218103 (2008)) came out 
on a related subject. These authors also consider the effect of fluctuations in the motion 
of bacteria. However, their goal is different. They are mainly interested in deriving transport 
coefficients from microscopic models, so they do not take into account the long-range interaction 
between bacteria due to chemotaxis. Alternatively, in our approach, the transport coefficients 

and x appearing in the Langevin equations are introduced phenomenologically but long- 
range interaction between bacteria due to chemotaxis is fully taken into account. Therefore, 
these two independent studies are complementary to each other. 
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